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Abstract 

The oscillation frequencies of collective excitations of a trapped Bose-Einstein 
condensate, when calculated in the mean-field approximation and in the 
Thomas- Fermi limit, are independent of the scattering length a. We calculate 
the leading corrections to the frequencies from quantum fluctuations around 
the mean field. The semiclassical correction is proportional to A rl / 5 a 6/ ' 5 , 
where N is the number of atoms in the condensate. The correction is positive 
semidefinite and is zero for surface modes whose eigenfunctions have a van- 
ishing laplacian. The shift in the frequency of the lowest quadrupole mode 
for an axially symmetric trap is large enough that it should be measurable in 
future experiments. 
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The achievement of Bose-Einstein condensation in trapped atomic vapors Jl| has opened 
up new possibilities for the experimental study of nonhomogeneous interacting Bose gases. 
This system can be conveniently described by a quantum field theory, but the quantum 
field equations are extremely difficult to solve in general. Fortunately, many of the basic 
properties of the condensates in existing experiments are described accurately by the mean- 
field approximation, which reduces the problem to solving the Gross-Pitaevskii equation, 
a classical equation for the mean field. These properties include the density profile of the 
ground state and collective oscillations of the condensate. Fluctuations of the quantum 
field around the mean field provide corrections to mean-field predictions that scale as the 
square root of the peak number density of the atoms. These corrections should be observable 
through precise measurements of the properties of condensates containing a large number 
of atoms. 

Among the most beautiful experiments on Bose-Einstein condensates are the excitation 
and study of collective oscillations of the condensates. Among the properties of the collective 
modes that have been measured are the frequencies of the lowest normal modes at extremely 
low temperatures , the temperature dependence of the frequencies and their damping 

rates, []5],[| , and the speed of sound || . The frequencies of the normal modes can be measured 
very accurately. In one experiment, they have been measured to better than a fraction of a 
percent ||. These frequencies are therefore an ideal place to look for the effects of quantum 
field fluctuations. 

Most previous studies of the collective excitations of a trapped condensate at zero tem- 
perature have been carried out using the mean-field approximation or the Bogoliubov ap- 
proximation, which is equivalent at T = 0. The collective excitations have been studied 
by variational methods [J], by numerical methods ||, and by analytic methods PJTOf. In 
the Thomas- Fermi limit, the oscillation frequencies for a condensate in a harmonic potential 
can be obtained analytically and are independent of the interaction strength of the atoms 
and the number of atoms in the trap |5j. Quantum fluctuations around the mean field give 
fractional corrections to the frequencies that are proportional to \fpcfi-, where p is the peak 
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number density and a is the S-wave scattering length of the atoms. The purpose of this 
paper is to calculate these corrections. 

A Bose-Einstein condensate containing a large number N of identical atoms in a trapping 
potential V(r) can be described by a quantum field theory with a single complex-valued field 
ip{r,t) that satisfies 

h 2 



ihip 



2m 



(1) 



where g = 16na. The chemical potential \x in (|I|) is to be adjusted so that the expectation 
value of J d 3 rip^ip, the total number of atoms, is equal to N. A collective oscillation of 
the Bose-Einstein condensate is described by a time-dependent state in the quantum field 
theory. Letting (. . .) denote the expectation value in that state, we can write the quantum 
field as ip = <p + ip, where <ft — ("0} is the mean field and ip is the quantum fluctuation field 
that satisfies (ip) = 0. The number density p and the current density j in that state are 



p=\<P\ 2 + 

j = -i-^ (0*V0 - 0V0*) - i^(ftV#-Vftj>). 



(2) 
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The continuity equation p + V • j = follows from the U(l) symmetry of the field theory. 
Taking the expectation value of ([[]), we obtain a partial differential equation for <fi that 
involves the matrix elements of ip: 
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In the mean-field approximation, the matrix elements are neglected and (U) reduces to 
the Gross-Pitaevskii equation. The effects of quantum field fluctuations are contained in 
the matrix elements, which can be expressed as functionals of the mean field (p. They are 
dominated by modes of ip with wavelengths comparable to the local coherence length £, which 
is the scale of the wavelength above which the excitations of the condensate are modified 
by collective effects. In the Thomas- Fermi limit, £ is much smaller than the length scale for 
significant variations in 0. The matrix elements can therefore be expanded in derivatives 



of by generalizing the methods developed in Refs. JTTjJT^] to the case of a time-dependent 
state. Including terms through first order in g 3 ^ 2 \(j)\ and through first order in derivatives of 
cf), we obtain |p~3f 
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To verify the continuity equation p + V • j = to first order in g 3//2 |0| and to second order in 
gradients of 0, we must include terms in (i/j^i/j) and (ipip) that are second order in gradients 
of 0. These terms are rather complicated |12| , but fortunately they are not needed to derive 
the equation for the time evolution of p in the Thomas- Fermi limit. By differentiating the 
continuity equation with respect to time, using (Q) to eliminate from the time derivative of 
j, and using @ to eliminate in favor of p, we can obtain a closed equation for p. Keeping 
only terms through first order in ((? 3 p) 1 ^ 2 and through second order in gradients of p, this 
equation is 
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This equation gives the time- dependence of the density in the Thomas-Fermi limit and 
to first order in the semiclassical expansion. The ground state number density po i n this 
approximation satisfies 



1 + 
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To deduce the equations for small-amplitude oscillations of the condensate, we substitute 
P = Po + P into © and linearize in p: 
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In the mean-field approximation, the ((? 3 po) 1 ^ 2 corrections in ( |TUD and ( [TT] ) are neglected. 
The solution for the ground state number density po vanishes at points r where V(r) > p 



and is given by po = 4m(p — V)/fi 2 g at points where V(r) < fj,. Assuming that p has the 
time dependence exp(iu)t), we obtain the eigenvalue equation H p = Cu 2 p, where 

H = -^V.(poV). (12) 

The "hamiltonian" Hq, whose eigenvalues give the mean- field approximations to the oscil- 
lation frequencies in the Thomas- Fermi limit, was first derived by Stringari |J. 

To obtain the semiclassical correction to the frequencies, we expand ( |TT| ) to first order 
in {g 3 po) 1 ^ 2 - The solution of ([K]) to first order in (g 3 p ) 1//2 is p = p [l — (g 3 ^) 1 / 2 /67r 2 }. 
Inserting this into ( |TT|) and assuming that p has the time-dependence exp(iut), we obtain 
an eigenvalue equation for u>. Absorbing a factor of 1 + g 3 ^ 2 PQ 2 /{^ 2 ) into p, our equation 
becomes the eigenvalue equation for a hermitian operator: (H + Hi)p = u 2 p, where 



ff ' = -wK v + v ^ /2 )- < 13 > 



The eigenvalues of H +Hi give the oscillation frequencies in the Thomas- Fermi limit through 
first order in the semiclassical expansion. The semiclassical corrections to the frequencies 
can be obtained by using the standard results of time-independent perturbation theory. If u 2 
is a nondegenerate eigenvalue of H with eigenfunction p n , then the semiclassical correction 



= uj2 n -u 2 n is 



ft 2 ^/ 2 Re/cfrpf p- ;V 2 p~ r 
487T 2 m 2 J d 3 r \p T 



6U n = --^T— 2 rM _^ ,2 • ( 14 ) 



If the eigenvalue u 2 is degenerate, then the expression (|14]) holds only in the basis consisting 
of states that diagonalize the restriction of Hi to the subspace of states with eigenvalue a) 2 . 
Note that the correction (|T4] ) is positive or zero. This follows from the fact that H\ is pro- 

—3/2 

portional to the anticommutator of the positive operator p and the positive-semidefinite 
operator —V 2 . The correction (|n|) is zero only if the laplacian of the eigenfunction p n 
vanishes. 

We first consider a condensate trapped in the spherically symmetric parabolic poten- 
tial V(r) = muo^r 2 /2. The ground state number density po in the mean-field approxi- 



mation and in the Thomas- Fermi limit vanishes outside a sphere of radius R given by 
N = m 2 uj 2 R 5 1 \lhh 2 a). Inside that sphere, it has the form 

Mr) = 2 -^A(R 2 -r% r<R. (15) 
h g 

The normal modes for oscillation of the condensate around the ground state can be la- 
belled by a radial quantum number n and angular momentum quantum numbers I and m. 
The mean-field frequencies u n e depend on the quantum numbers n and £, and were first 
determined by Stringari ||: 
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2n 2 + 3n + (2n + 1)1 to 2 . (16) 



The n = modes are surface modes and their eigenf unctions are r e Yi m (8, 0). Since they 
have a vanishing laplacian, the semiclassical correction to the frequencies of these surface 
modes vanishes: 5u}q £ = 0. The eigenfunctions for n = 1 are (R 2 — j^r 2 )r e Yi m (9, (j>) for 



r < R. Inserting this into (H|), we find that the shifts in the frequencies are 



where (z) n = T(z + n)/T(z) is the Pochhammer symbol. The fractional shift in the frequency 
is proportional to muoaR/h, which scales like iV 1 It therefore increases very slowly 

with the number of atoms in the condensate. 

We next consider a condensate trapped in the axially symmetric parabolic potential 
V(p,z) = m(uj'j_p 2 + uo 2 z 2 )/2. The ground state number density in the mean-field approxi- 
mation and in the Thomas-Fermi limit takes the form 

Po(r) = 2 -^(R 2 - p 2 - AV), p 2 + \ 2 z 2 < R 2 , (18) 
k g 

where A = uh/ux and R is given by iV = m 2 uj 2 _R!' /(15h 2 a\). The eigenfunctions have 
the form p\ m \ e im<i > P {z , p 2 ), where m is the azimuthal angular momentum quantum number 
and P is a polynomial in z and p 2 . The number of independent solutions for which P 
is a k'th order polynomial in z is int(l + k/2). For k = and k — 1, the frequencies 
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in the mean-field approximation are u 2 m = \m\uj_ and uf m = (\m\ + \ 2 )ujj_, respectively. 
The corresponding eigenfunctions are p\ m \e. im <f> and 

zp \m\ e imcf>_ gince 

these eigenfunctions 

have vanishing laplacian, the semiclassical corrections to the frequencies vanish: Su 2 m = 
Suf m = 0. For each value of m, there are two independent eigenfunctions with k = 2. Their 
frequencies in the mean-field approximation are 



-2 



2m± 



(2{\m\ + 1) + §A 2 ± v /(|m|+2) 2 -(|m|+4)A 2 + |A 4 ) u\. (19) 



These were first derived by Stringari for m = M and by Ohberg et al. and by Fliesser et 
al. for general m JTD] . The eigenfunctions have the form (R 2 + a m ±p 2 + p m ±\ 2 z 2 )p^e vm ^, 



where a m ± and j3 m ± are complicated functions of \m\ and A. Using (0), we find that the 

semiclassical corrections to the frequencies are 

2 3(i)| m | +4 /±(H, A) fmuj ± aR\ 2 

^ = V2(\m\ + 4)l 

u l , , 1 , lx2 , (|m| + l)(H+2)-iA 2 + |A 4 
/±( m L A) = m + 1 + ^A ± . =. (20) 

v /(|m|+2) 2 -(|m|+4)A 2 + |A 4 
If we set A = 1, this reduces to ( |TTD with t — > |m|. 

The most accurate measurement to date of an oscillation frequency for a Bose-Einstein 
condensate was carried out by Stamper-Kurn et al. |||. They trapped iV = 1.5 X 10 7 atoms 
of 23 Na, whose scattering length is a = 2.75 nm, in an axially symmetric potential. The axial 
trapping frequency was measured to be cjii/27r = 16.93(2) Hz and the transverse trapping 
frequency was estimated to be u^/2n = 230 Hz, so that A ~ 0.074. They measured the 
ratio W2o-/^|| °f the lowest quadrupole oscillation frequency to the axial trap frequency 
to be 1.569(4). The fractional uncertainty in the measurement is 0.25%. The mean-field 
prediction from (|19|) is u^a-I^W — 1-5806, which is larger by about 3 standard deviations. 



Including the semiclassical correction (20), we obtain the prediction 1.5813. The fractional 



increase is 0.04%, which is about g of a standard deviation. Thus there is a several standard 
deviation discrepancy between the prediction and the measurement. 

In order to assess the significance of the discrepancy, we must consider other sources of 
corrections to the frequency, including edge effects, nonlinear effects, and nonzero temper- 



ature effects. The correction from edge effects arises from the breakdown of the Thomas- 
Fermi approximation near the edge of the condensate. The fractional correction scales like 
(h/mooR 2 ) 2 which is approximately 4% for the experiment in Ref. ||, so these cor- 



rections are probably not negligible. This source of theoretical error can be eliminated by 
avoiding the Thomas- Fermi limit, and instead calculating the frequencies numerically by 
solving the mean-field equations. The frequency shift from nonlinear effects has been stud- 



ied by Dalfovo et al. ||15|| . The fractional shift for the lowest quadrupole mode was found to 
be +0.09A 2 , where A is the relative magnitude of the oscillation in the length of the con- 
densate. Taking the value A « 0.1 from Ref. ||, we find that the fractional frequency shift 
is 0.9%. This source of error can be eliminated experimentally by decreasing the amplitude 
of the excitation of the condensate. It was suggested in Ref. || that the discrepancy in the 
frequency can be explained by nonzero temperature corrections, which have the correct sign. 



While there have been studies of the temperature dependence of the frequencies ||16|| , there 
don't seem to be any quantitative treatments for the low temperatures appropriate to the 
experiment of Ref. || . 

We have calculated the semiclassical corrections to the frequencies for small oscillations 
of a Bose-Einstein condensate in the Thomas- Fermi limit. The frequency shifts are less than 
an order of magnitude smaller than the uncertainties in some previous frequency measure- 
ments. These shifts should therefore be measurable by more accurate experiments. Accurate 
theoretical predictions require going beyond the Thomas- Fermi limit in calculating the semi- 
classical frequency shifts. Separating the semiclassical correction from that due to nonzero 
temperature may require the development of new techniques for measuring the tempera- 
tures of extremely cold condensates. A measurement of the semiclassical frequency shift 
in agreement with theoretical predictions would mark a new threshold in our quantitative 
understanding of Bose-Einstein condensed gases. 

While this work was being prepared for publication, a preprint by Pitaevskii and Stringari 
T7| appeared in which the authors also calculate the semiclassical corrections to the os- 



cillation frequencies of a trapped Bose-Einstein condensate. Their results are in perfect 



agreement with ours. 

We thank J.O. Andersen for valuable discussions. This work was supported in part by the 
U.S. Department of Energy Division of High Energy Physics and by a Research Experience 
for Undergraduates site grant at the Ohio State University. 
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